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We propose the model for magnetic properties of the magnetic super- 
conductor RuSr2GdCu208, which incorporates the theory of the supercon- 
ducting/ferromagnetic multilayers. The transition line Tfi{h), on which the 
Josephson coupled superconducting planes are decoupled, i.e. jc{Td) = 0, 
is calculated as a function of the exchange energy h. As the result of this 
decoupling a nonmonotonic behavior of magnetic properties, like the lower 
critical field Hd, Josephson plasma frequency, etc. is realized near (or by 
crossing) the Td{h) line. The obtained results are used in analyzing the 
newly discovered antiferromagnetic ruthenocuprate RuSr2GdCu208 with 
possible weak ferromagnetic order in the RuO planes. 
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I. INTRODUCTION 



The physics of magnetic superconductors is interesting due to competition of magnetic 
order and singlet superconductivity in bulk materials. The problem of their coexistence 
was first set up theoretically in the pioneering work by V. L. Ginzburg in 1956, while 
the experimental progress in the field begun after the discovery of ternary rare earth 
(RE) compounds (RE)Rh4B4 and (RE)Mo6X8 (X=S,Se) with a regular distribution of 
localized RE magnetic moments. It turned out that in many of these systems supercon- 
ductivity (with the critical temperature Tc) coexists rather easily with antiferromagnetic 
(AF) order (with the critical Neel temperature Tn), where usually the situation with 
Tjv < Tc is realized . Due to their antagonistic characters singlet superconductivity and 
ferromagnetic order cannot coexist in bulk samples with realistic physical parameters. 
However, under certain conditions the ferromagnetic order is transformed, in the presence 
of superconductivity, into a spiral or domain-like structure - depending on the type and 
strength of magnetic anisotropy in the system As the result of this competition, these 
two orderings coexist in a limited temperature interval Tc2 < T < Tm (the reentrant 
behavior) in ErRh4B4 and HoMoeSg, or even down to T = OK in HoMoeSeg, where Tm 
is the critical temperature for the existence of the inhomogeneous magnetic order. The 
coexistence region in ErRh4B4 is narrow where Tc = 8.7 K, Tm ~ 0.8 K, Tc2 ~ 0.7 K, while 
for HoMoeSs it is even narrower with Tc = 1.8 K, Tm ~ 0.74 K, Tc2 ~ 0.7 K - see Refs. 
0,0. In most of the new quaternary rare-earth compounds (RE)Ni2B2C the antifer- 
romagnetic order and superconductivity coexist up to T = OK |@, while in HoNi2B2C 
an additional oscillatory magnetic structure is realized in a limited temperature interval. 
This oscillatory magnetic structure competes strongly with superconductivity giving rise 
to reentrant behavior in this compound |^. Recently the Pobell's group in Bayreuth 
made a remarkable discovery of the coexistence of superconductivity and nuclear mag- 
netic order in Auln2 with Tc = 0, 207 K and Tm = 35/iK. This exciting phenomenon was 
explained in Ref. where it is argued that superconductivity can coexist either with 
spiral or domain-like nuclear magnetic ordering only, depending on the strength of mag- 
netic anisotropy in this cubic system. Important contribution to the physics of magnetic 
superconductors has been done in Ref. 0, where for the first time was proposed the co- 
existence of weak-ferromagnetism and superconductivity. In such a case the spontaneous 
vortex state due to weak-ferromagnetism is also possible. 

We point out that in the above cited magnetic superconductors the exchange inter- 
action (between localized magnetic moments and conduction electrons) influences super- 
conductivity much stronger than the electromagnetic interaction. The latter is due to the 
localized magnetic moments which create dipolar magnetic field, thus affecting the orbital 
motion of superconducting electrons. 

Recently, a new class of magnetic superconductors based on layered perovskite 
ruthenocuprate compound RuSr2GdCu208 comprising Cu02 bilayers and RuO mono- 
layers has been synthesized This compound belongs also to the class of high-Tc 
superconductors (HTS). A subsequent study of transport and magnetic properties has 
revealed that it exhibits some kind of ferromagnetic order at the critical temperature 
T/v = (133 — 137) K. The polarized neutron scattering measurements [|10| show that the 
magnetic structure (which appears at T/v) is predominantly antiferromagnetic with a Ru 
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magnetic moment /xru ~ l.lSfiB along the c-axis at low temperature. The same mea- 
surements put an upper limit ~ O.lfiB to any net ferromagnetic zero-field Ru moment. 
Concerning the last point the important results came from magnetization measurements 
first reported in Ref . , which show a hysteresis loop and remanent magnetization. The 
latter hints to existence of a ferromagnetic component in the system. Recent magne- 



tization measurements on RuSr2EuCu2 08 |TT]] give evidence for a small ferromagnetic 



component, which lies probably parallel to the RuO plane, with the magnetic moment 
(per Ru) ~ O.OSyUB at 5K consistent with the neutron scattering data [TU|. Note that 
the smaller value of magnetic moment (O.OSyUs) in this compound tells us that in the Gd- 
compound some admixture of the large Gd moment might take place. This conclusion is 
also confirmed by the zero-field muon spin rotation (ZF-/iSR) measurements |]12[ which 
provide important evidence that the magnetic order is homogeneous on a microscopic 
scale and accounts for most of the sample volume. At lower temperatures the supercon- 
ductivity sets in at Tc = (35 — 45) K without affecting the AF order |jTO|, [l^ notably. 
This fact means that superconductivity - which is realized predominantly in the Cu02 
planes, and magnetic order - which is present only in the RuO planes, interact rather 
weakly, i.e. these two orders are separated spatially. Recently, it was reported that 
in Rui_a;Sr2GdCu2+a;08-j, the highest superconducting critical temperature reaches 72 K 
for X = 0.3 — 0.4, while there is no sign of the weak-ferromagnetic (WF) component in 
the RuO planes. 

It seems that RuSr2GdCu208 has very interesting magnetic properties, which might 
result in the absence of Meissner phase in some samples [0] , , while in some others it 
is realized 0, [jl6|. (This problem will be briefiy discussed in Section IV.) 

In this paper we propose a model of layered magnetic superconductor with weak- 
ferromagnetism, which might be relevant for the RuSr2GdCu208 compound - the SWF 
model. This model, studied in Section II, assumes the existence of S/F multilayers with 
small hopping parameter t between S (superconducting) and F (ferromagnetic) planes 
along the c-axis , i.e. t < Tc. As a result the small t gives rise to an effective Josephson 
coupling current jc between superconducting planes. It turns out that jc is suppressed by 
the exchange field - present in the F-plane only, which causes drastic changes of magnetic 
properties. The Gibbs free-energy Q of such a magnetic superconductor with both AF 
and WF orderings in external magnetic field H is formulated in Section III. Based on it 
the lower critical field Hd is also studied there. The estimation of theoretical parameters 
of the SWF model from the experimental results in RuSr2GdCu208 is done in Section IV, 
where the obtained results are discussed too. 



II. MODEL FOR S/F ATOMIC MULTILAYER AND JOSEPHSON CURRENT 

As was mentioned above we consider the magnetic superconductor RuSr2GdCu208 as 
a prototype for S/F atomic multilayers by assuming good conduction in Cu02 planes - 
with the quasiparticle spectrum ^s{p), and a small hopping parameter t <Tc between the 
S- and F-planes (i.e., along the c-axis). The second assumption is related to the existence 
of WF order (with the magnetization M lying in the RuO planes) which gives rise to 
an effective exchange field parameter h = heat- The latter affects spins of conduction 
electrons with the dispersion ^f(p) moving in the normal conducting RuO planes. (The 
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afe-plane is sometimes labeled by the xy-plane.) The parameter h can be related to an 
effective spontaneous spin S^s (magnetization normalized to saturation magnetization) in 
the a6-plane, i.e. h = J'^'^Ses - see also Section IV. 

The electronic part of the SWF model is similar to the model in the Ref. and in 



what follows the same notation is used. According to this model the elementary cell of the 
superlattice consists of one superconducting and one ferromagnetic layer which are both 
metallic. For simplicity it is supposed here that both layers have similar quasiparticle 
energy spectra, i.e. ,^(p) = (^s(p) ~ ^f(p))- It is also assumed, that the superconduc- 
tivity is realized in S-planes (Cu02 planes) with pairing coupling g{p) (having in mind 
apphcation to the HTS compound RuSr2GdCu208 the clean limit, I, is supposed). 
The Hamiltonian of the system is given by 

H = Hq + Hintl + Hint2i 

Ho= ^(P)^U,M^n,i,aiP) + 4,l,<x(P)«n,-l,<x(p) + 4+l,-l,a(p)«n,l,<x(p) + h.C. 

p,n,i,a 



Hintl = \ Yl 3^^^ ~ P2)4,l,<x(Pl)4,l,-<x(-Pl)«n,l,-a(-P2)a„,i,^(p2) 
pi,P2,n,o- 

Hint2 = ^^4,-l,a(p)V-l,'x(P)) 



p,n,cr 



where ajjjo-(p) is the creation operator of an electron with spin a (the quantization axis 
is parallel to the a6-plane) in the n-th elementary cell and momentum p in the layer i is 
parallel to the a6-plane, where i = 1 for the S layer, and i = —1 for the F layer. Since the 
obtained results below are qualitatively similar for s- and d-wave pairing, the calculations 
were done for s-wave pairing where ^^(p) = go is constant, while quantitative changes due 
to d-wave pairing are discussed below and in Section IV. 

By assuming that the order parameter changes from cell to cell in the manner A„ = 
I Ale**^" (with = kn in absence of orbital effects) the quasiparticle Green's functions are 



obtained in the standard way . The self-consistency equation for the order parameter 
|A| reads [O] 



where A = gop{0) and p(0) = m\\/2Tx is the electron density of states at the Fermi level in 

the normal state, uj± = iu ± ^{p)i ^± = uj± + uj = 7cT{2n + 1). The quasimomentum 

q lies in the direction perpendicular to the layers, and Tg = 2t cos(g/2)e*''/^ and 7^+^ = 
2t cos((g + fc)/2)e*(«+'^)/2_ 

The free energy JF in the superconducting state is obtained by using the following 
relation 

where the expression for F^i is obtained in [17 
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In order to study transport and magnetic properties in magnetic field we need to know 
the supercurrent flowing across the layers (along the c-axis in RuSr2GdCu208). In this 
case, the vector potential = AzGz enters the Hamiltonian through the substitution t 
^^±tedA^/c^ where d is S-F interlayer distance, and the part of the Hamiltonian depending 
on Az is given by 

Ha=Y. i[<iM<^n,-i,Me'''^'''' + ai^i,-iAvWxAv)e'''^''' + h.c]. (3) 

p,n,i,o- 

The supercurrent across the planes is obtained by the standard procedure 

2d 5Az ■ ^ ' 

Note that the Josephson supercurrent in the S/F superlattice is carried by Andreev bound 
states, similarly to the S/N case. In S layer the supercurrent is carried by Cooper pairs, 
but in N layer it flows via quasiparticles, which recondense in the next S layer; bound 
states represent this process |jl8[ . 



In the case of small hopping parameter t <^Tc the Josephson current along the c-axis 
is obtained in leading order (proportional to t^) by standard perturbation theory. After 
the integration over the energy ^ it reads 

,. = 4e.|APiV(0)TX:{2. J^fj^^J:^) (5) 



= 3c sin k, 

where R{ijj) = (| Ap + /i^)^ + 4/i^cc'^ and Q{uj) = y^u'^ + |A|2. As in the standard Josephson 
effect, the supercurrent jz is proportional to sin/c, k being the phase difference between 
n- and + 1 S- layers. 

In what follows we calculate numerically the critical current jc (in Eq. (^) ) at any point 
of the phase diagram (T, h) by replacing |A| — Aq(T), where Ao(T) is given by the BCS 
theory. The latter is correct due to the smallness of t, in which case Tc is practically 
unaffected by the exchange fleld, i.e. Tc ~ Tco up to the second order terms in t/TcQ. 

From Eq.(|]) it comes out in particular, that near and for h = one has jc > 0, while 
jc < for h ^ Tc. The change of sign oi jc (near Tc), which corresponds to the transition 
from = to = TT in the ground state, occurs at he = 3.77Tco, in accordance with the 
calculation in Jl^. At low temperatures, T — > 0, jc goes to zero at /z/Ao(0) ^ 1/2, which 



just corresponds to hco = 0.87Tco at T = 0, again in accordance with Ref. IjT^. Note 
that the same approach if applied to d-wave pairing [1^ gives = O.QTco at T = 0. 



The sign-change of jc is related to the transition from the "0" -phase to "vr" -phase. This 
transition goes smoothly if we take into account the higher order term (~ cos 2k) in the 
free-energy, in fact it means that the width of the region Ah where the transition from 
"0" -phase to "vr" -phase occurs is of the order of Ah ~ t^/T^Q- In the case of weak hopping 
t <^ Tco this region is very narrow and we may deflne the decoupling line jc{Td, h) = 
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which results in the (T^, h) phase diagram shown in Fig.l. 




FIGURE l.The (T, /i)-phase diagram for the case t <C Tco- The Josepson decouphng hne jc{T, h) = - 
black hne 

The temperature dependence of the Josephson penetration depth Aj = 
a/c^o/Stt^ \jc\ {2d) [^1, where 0o is the flux quantum, is shown in Fig.2 for various h 0. 
Here 2d is the period of the multilayer. One should note its nonmonotonic behavior if 
h ^ 0, particularly when h ~ T^. Based on these results one can analyze some magnetic 
properties, like the lower critical field Hd in the a6-plane. 
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FIGURE 2. The T-dependence of the Josephson penetration depth Aj(T, h) for various h. We defined 
Ao = UeTcO, jcQ = ep{0)t*/Al, and Ajo = ^/c(f>o/^6^^^djco■ 



III. GIBBS ENERGY AND IN-PLANE CRITICAL FIELD Hci 
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A. Gibbs energy 



In order to calculate the lower critical field Hd (and the possible absence of Meissner 
phase |T^, |]T3|) we need the Gibbs energy functional Q. Having in mind the application to 



the RuSr2GdCu208 we assume, according to the neutron scattering data [T^, that in the 
magnetic subsystem (F-layers coinciding with RuO planes) AF order with spins along the 
c-axis is realized at Tjy ^ T^. The AF order parameter is L = L^e^. The magnetization 
measurements 0, [|ll|] imply WF order with the magnetization lying (most probably) 
parallel to the a6-planes and with the effective moment, |M| /nRu = fics < O.lfis, where 
riRu is the density of the magnetic Ru ions. In order to construct the magnetic free- 
energy it is necessary to know the symmetry of the system as well as orientations of the 
easy axis in different sublattices. So, for instance if the magnetic anisotropy energy on 
different sublattices are unequal, then one expects the WF order to be realized. However, 
at present there are no sufficient experimental data on the local lattice distortion which 
might favor WF order and accordingly the preferred direction of M. The above discussed 
experiments [|10|, ||12|, |0, ||1T| suggest only that M is in the a6-plane, i.e. M = M^e^, 
where is in the a6-plane. In the following we define the x-axis by = e^. As the 
result of this analysis the SWF model contains the following order parameters: Lz for the 
AF order, for the WF order and An{x,y) for the S order. 

In the applied magnetic field H the Gibbs energy of the layered magnetic supercon- 
ductor reads - see also 0,0], 

g [A„, L, M, B; H] = 1 dV (^J^m [L, M] + - ^^M)' _ ^ 

+ ^(2rf) f dxdyJ^s[^n,A], 

where A„ = A„(x, y), and L, M, B are also coordinate dependent. The magnetic field 
B = rot A is due to the dipolar field created by the magnetic moments, the external 
magnetic field and the superconducting screening current. The vector potential A = 
Aab+Ac = Aab^ab + A^^z coutaius the component, Aab, in the a6-plane, and Az along the 
c-axis. 

The magnetic free-energy density functional J^m [L, M] which mimics the experimental 
results in RuSr2GdCu208 with L = L^e^ and M = M^e^ is given by the following 
phenomenological expression 



Tm [L, M] = + ^L^ + ^M^ - ^LzM, + F„(L, M) + ^^(VL, VM) + ... 



(7) 



Although this expression is quantitatively correct near the AF transition at T/v it is also 
suitable for semiquantitative analysis even below superconducting transition temperature 
Tc, due to the smallness of M and 7. The first two terms describe the AF order (a = 
a'{T — T/v) < 0, /3 > 0), the third (5 > 0) and fourth (~ 7) terms describe the induced 
WF by the AF order. The parameters a and 5 are due to the exchange interaction 
(between Ru spins in RuSr2GdCu208), where one has a' ~ l/Oem and 5 ~ T^jdem with 
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Oem = 27r/i^ ~ 1 K- The unknown anisotropy term J^a fixes the direction of L and M, i.e. 

Since in the following we analyze the lower critical field along the a6-plane, with 
characteristic length-scales Aa6,Aj 3> ^afe, c?Ru-Ru, where ^ab is the coherence length and 
c^Ru-Ru is the Ru-Ru distance, it is justified to omit the gradient term jFg(VL, VM). By 
minimizing [L, M] with respect to and one gets (at temperatures < T < T^) 

Ml = (8) 

The neutron scattering and magnetization measurements give limits for (L^/riRu) ~ (1 — 
1.2)//^, and (M°/nRu) < (0.05 — O.l)/!^, which implies an upper limit for 7, i.e. (7/^) < 
(0.05-0.1). 

According to experiments |12| the AF (and WF) ordering is practically unaf- 
fected by the appearance of superconductivity, then it is reasonable to neglect the effect 
of superconductivity on the exchange (RKKY) interaction in jFg [A„, A]. Therefore we 
keep in JF5 [A„, A] the electromagnetic interaction between superconducting electrons and 
magnetic order only 

:Fs [A„, A] = [|A„| , 0] + (^^y + 1^(1 - cosx„,„+i), (9) 

where Fs'[|A„|,0] is the condensation energy and A„ = | A„| exp («(/)„). Note that the 
exchange interaction between conduction electrons and localized Ru moments affects su- 
perconductivity by renormalizing jc, which is function of h. 
The current in the a&-plane jah reads 

iab = --Arri^ab - P-^ab^n), (10) 

where Xab is the bulk London penetration depth in the ab superconducting layers (we 
assumed = A;, = Xab)- The last term depends on the gauge invariant phase Xn,n+i 

Xn,n+1 = - 7 , (11) 

00 

which characterizes the effective Josephson coupling between two neighboring S-planes 
with the distance 2c?. It is due to the hopping between S- and F-planes and jc = jciT, h) 
is determined by Eq.(H). 



B. Lower critical field Hd 

Let us calculate the lower critical field for the case when the magnetic field H and 
the single vortex are along the magnetization M = M^Gx-, i-e. H = H^Gx and B = Bxe.x- 
By the standard minimization procedure of the Gibbs free-energy Q [A„, L, M, B; H] with 
respect to A„, L, M, B, and by assuming the continuum limit |]2D|, one gets the complete 
set of equations for these quantities as well as the Gibbs free-energy of the vortex - see 
also i, i. 
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(5 + 47r)M,-7L,-5,. = 0. 



(12) 



The Maxwell equation for the magnetic field B reads 

4:71 

rot(B-47rM) = — j„ (13) 

where 

h=iab+iz- (14) 

The in-plane current jab is given by Eq. ([ToD while is the Josephson current between 
planes 

iz = jcezsmxn,n+i- (15) 

The phase Xn,n+i is given by Eq.([TTD. In the following we assume that the vortex axis, 
B and the external field H are along the x-axis. By the standard procedure we get the 
equation for the single vortex (centered on the origin) 

The parameter p = 6/{6 + An) takes into account the additional screening due to the 
appearance of the WF order. After straightforward transformations the Gibbs energy of 
the vortex, Q^, has the form 

where 

M° is approximately given by Eq.(H). 

The Eq.(|T6[) can be generalized by taking into account the nonlinear core effects | pO|| . 
In that case the solution is 

BAy,z)= ^ \' ■ Koi-). (18) 



where R = \/2/^/Aj + z'^/X^f^ and Kq is the Bessel function of the zeroth order of an 
imaginary argument. Inserting such a solution into Eq.(|l7|) a straightforward calculation 
gives the lower critical field Hd from the condition Qv = 

H..t + A^pMl = Hd ^ In (19) 

where M° = {■y/6)Lz. We stress that the logarithmic factor in Eq.(|TUp is due to the 
nonhnear core effects of the Josephson vortex pO| . Note that in systems with T/v ^ 9em , 
like in RuSr2GdCu208, one has p ~ 1. 
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From (1T9|) it is seen that for 



^0 



In 



Xab^/P 



d 



(20) 



spontaneous vortices appear in the system. This condition is more easily reahzed 
near the "0"-to-"7r" transition (decouphng) hne Td{h), i.e. when Aj is significantly in- 
creased. This means that in systems where the exchange parameter fulfills the condition 
0.87Tco < h < 3.77Tco for s-wave pairing, while for d-wave pairing 0.87Tco is replaced by 
0.6Tco, then by lowering the temperature the Hci{T, h) shows pronounced nonmonotonic 
behavior reaching minimum at the "0"-"7r" boundary line as it is seen in Fig. 3. 
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FIGURE 3. The T-dependence of the lower critical field, Hci(T,h), in the a6-plane for various h.We 
defined Hcw = (p$o/47rAjoAafc) In {Xab^/p/d) . 



IV. COMPARISON WITH THE EXPERIMENT AND DISCUSSION 

Let us discuss some relevant points related to the interpretation of the obtained results 
on RuSr2GdCu208. 

(i) In order to analyze the magnetic properties the value of the exchange field parame- 
ter h = J°'^Sf.fi is needed. If one takes the experimental value S^s ~ 0.1 one gets 7/5 ~ 0.1 
since Ses ~ 7/^- However, at present we do not know the relation between J"'^ and T/y. As 
^Ru-Ru '^Ru-Ru expects that the coupling of spins along the c-direction, J'^ is much 
smaller than along the a6-plane, J""^ . In such a situation one has T/v ~ J"^/ ln( J"^/ J^). 
Then, J""^ > T/v and /i ~ 10 — 20K. We pay attention that there are evidences that in 



the underdoped HTS materials d-wave pairing is realized |21|. In that case the point on 
the phase diagram jc{T = 0, he) = is realized for h^^^ = 0.6Tc. According to the specific 
heat measurements [0 in RuSr2GdCu208 with = 30 — 40K this compound behaves 
like an underdoped HTS material. If it is so it gives h^cf^ ~ 20 K, i.e. h is near to he and 
a nonmonotonic behavior of Hd is expected as shown in Fig. 3. 
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(ii) If M° fulfills Eq.(pO|) then there is a spontaneous vortex state and the Meissner 
effect is absent. In opposite case the Meissner state is realized. 

(iii) It may happen that < Hd in some temperature intervals and > Hd in 
the interval between, which case corresponds to a reentrant behavior. 

(iv) At present the origin and the magnitude of the parameter 7 in Eqs.(^,^ is un- 
known. However, it may also happen that in polycrystalline samples strains induce addi- 
tional changes of this quantity. A drastic case might be realized if the symmetry of the 
crystal implies that 7 = 0. Even in that case strains in samples, for instance the compo- 
nent axy, can induce a magnetic moment in piezomagnetic systems, i.e. M° ~ CxyLz thus 
producing weak ferromagnetism. If strains in a sample are such that M° > Hd then the 
Meissner phase is not realized as reported in In such a way one could reconcile 



the opposite claims on existence , and nonexistence [jT^ , [jT5| , of the Meissner phase 
in differently prepared samples of RuSr2GdCu208. 

(v) Based on the above analysis one expects that dynamical properties of such a system 
are very exotic. For systems near the decoupling line jc{Td, h) = there is a significant 
reduction of the Josephson plasma frequency ujq^sf ~ Vie ^o,jj I^Ml (compared to 
standard Josephson junction with tu^jj) for the waves propagating along the xy-planes 
in the S/F superlattice 

2 2 22 

^SF = ^0,SF + "^SfQ ■ 

Due to the microscopic character of the S/F superlattice one expects that Vgp ^ Vjj where 
vjj is the phase velocity for the Josephson junction made from bulk superconductors. This 
means that in a S/F superlattice, like for instance in RuSr2GdCu208, it is possible to tune 
lVq gp nonmonotonically and also to extract the radiation with much higher intensity than 
in the single Josephson junction. This is a matter of future investigations. 

In conclusion, we have shown that in a S/F superlattice with the exchange field h ^ Tc 
acting in F-planes only a nontrivial and nonmonotonic behavior of magnetic properties, 
like the lower critical field Hd, is realized. This property is due to the decrease of the 
effective Josephson coupling between S-planes by increasing h. 
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